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Abstract 

Within the hght-front framework, form factors for P ^ P and P ^ V tran- 
sitions (P: pseudoscalar meson, V: vector meson) due to the valence-quark 
configuration are calculated directly in the entire physical range of momentum 
transfer. The behavior of the form factors in the infinite quark mass limit are 
examined to see if the requirements of heavy-quark symmetry are fulfilled. 
We find that the Bauer-Stech-Wirbel type of light-front wave function fails 
to give a correct normalization for the Isgur-Wise function at zero recoil in 
P — > y transition. Some of the P ^ V form factors are found to depend on 
the recoiling direction of the daughter mesons relative to their parents. Thus, 
the inclusion of the non- valence contribution arising from quark-pair creation 
is mandatory in order to ensure that the physical form factors are independent 
of the recoiling direction. The main feature of the non-valence contribution is 
discussed. 

PACS numbers: 13.20, 14.40.J 
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I. INTRODUCTION 



The hadronic matrix element of weak P ^ P transition (P: pseudoscalar meson) is de- 
scribed by two form factors, whereas in general it requires four form factors to parametrize 
the weak matrix element ioi P ^ V transition {V: vector meson). Heavy quark symmetry 
predicts that, all the mesonic form factors in the infinite quark mass limit rriQ — > oo are re- 
lated to a single universal Isgur-Wise function [Q. The symmerty breaking l/mg corrections 
can be studied in a systematic framework, namely the heavy quark effective theory (for a 
review, see [^). The Isgur-Wise function is normalized to unity at zero recoil, but otherwise 
remains unknown. Phenomenologically, the hadronic form factors can be evaluated in various 
models among which the quark model is a popular one. However, since usual quark-model 
wave functions best resemble meson states in the rest frame or where the meson velocities 
are small, hence the form factors calculated in non-relativistic quark model or the MIT bag 
model are trustworthy only when the recoil momentum of the daughter meson relative to 
the parent meson is small. 

As the recoil momentum increases (corresponding to a decreasing g^), we have to start 
considering relativistic effects seriously. In particular, at the maximum recoil point = 
where the final meson could be highly relativistic, there is no reason to expect that the 
non-relativistic quark model is still applicable. A consistent treatment of the relativistic 
effects of the quark motion and spin in a bound state is a main issue of the relativistic 
quark model. To our knowledge, the light-front quark model [^,^ is the only relativistic 
quark model in which a consistent and fully relativistic treatment of quark spins and the 
center-of-mass motion can be carried out. This model has many advantages. For example, 
the light-front wave function is manifestly Lorentz invariant as it is expressed in terms of the 
momentum fraction variables (in "-|-" components) in analog to the parton distributions in 
the infinite momentum frame. Moreover, hadron spin can also be correctly constructed using 
the so-called Melosh rotation. The kinematic subgroup of the light-front formalism has the 
maximum number of interaction-free generators including the boost operator which describes 
the center-of-mass motion of the bound state (for a review of the light-front dynamics and 
light-front QCD, see §). 

The light-front quark model has been applied in the past to study the heavy-to-heavy 
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and heavy-to- light weak decays form factors However, the weak form factors were 

calculated only for < 0, whereas physical decays occur in the time-like region < < 
{Mi — MfY, with Mij being the initial and final meson masses. Hence extra assumptions 
are needed to extrapolate the form factors to cover the entire range of momentum transfer. 



In [|To| an ansatz for the dependence was made to extrapolate the form factors in the 
space-like region to the time-like region. Based on the dispersion formulation, form factors 
at > were obtained in |ll| by performing an analytic continuation from the space-like 
region. Finally, the weak form factors for P ^ P transition were calculated in 1^141 
for the first time for the entire range of g^, so that additional extrapolation assumptions 



are no longer required. This is based on the observation |15| that in the frame where the 
momentum transfer is purely longitudinal, i.e., q± = 0, g^ = g^g^ covers the entire range of 
momentum transfer. The price one has to pay is that, besides the conventional valence-quark 
contribution, one must also consider the non- valence configuration (or the so-called Z-graph) 
arising from quark-pair creation from the vacuum (see Fig. 1). The non- valence contribution 
vanishes if q~^ = 0, but is supposed to be important for heavy-to-light transition near zero 
recoil P,[l0|,|l5|,[l6| . Unfortunately, a reliable way of estimating the Z-graph contribution is 
still lacking. 

In the present paper we calculate the P ^ V form factors directly at time-like momentum 
transfers for the first time. We then study the mesonic form factors in the infinite quark mass 
limit to check if the light-front model calculations respect heavy quark symmetry. We are 
able to compute the Isgur-Wise function exactly since the non-valence contribution vanishes 
in the heavy-quark limit. It turns out that not all light-front wave functions give a correct 
normalization for the Isgur-Wise function at zero recoil in P —* V decay. In other words, the 
requirement of heavy quark symmetry can be utilized to rule out certain phenomenological 
wave functions. 

Another issue we would like to address in this work has to do with reference frame 
dependence of form factor. For a given g^, one can choose whether the recoiling daughter 
meson moves in the positive or negative z-direction relative to the parent meson, which we call 
the "+" and "— " reference frame, respectively. For some form factors in P ^ V transition, 
namely Aq, Ai, V, valence-quark and non- valence contributions are separately dependent 
on the choice of the "+" or "— " frame, but their sum should not. This demonstrates the 
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fact that it is mandatory to take into account the non- valence configuration in order to have 
physical predictions for form factors. This issue will be discussed in more details in Sections 
IID and IVC. 

This paper is organized as follows. In Section 11, the basic theoretical formalism is given 
and form factors for P ^ P and P ^ V transitions are derived. Section 111 is devoted to 
the discussion of the Isgur-Wise function. Numerical results are present and discussed in 
Section IV, and finally a summary is given in Section V. 



We will describe in this section the hght-front approach for the calculation of the weak 
mesonic form factors for pseudoscalar-to-pseudoscalar and pseudoscalar-to- vector transitions. 
The hadronic matrix elements will be evaluated at time-like momentum transfers, namely 

the physically accessible kinematic region < < Q'^ax- 

A meson bound state consisting of a quark qi and an antiquark q2 with total momentum 
P and spin S can be written as 



II. FRAMEWORK 



\MiP,S,S,)) = I {d'p,}{d'p2} 2{27i f5%P-p,-P2) 

X J2 ^^^^(pi,P2,Ai,A2) |gi(pi,Ai)g2(P2,A2)), 



(2.1) 



Ai,A2 



where pi and p2 are the on-mass-shell light-front momenta. 




(2.2) 



and 




k(pi,Ai)g(p2,A2)) = 6Upi)4.(P2)|0), 

{by{p'), b{{p)} = {dy(p'), diip)} = 2(27r)^ 5'{p' - p) 5vA. 



(2.3) 



In terms of the light-front relative momentum variables (x, k±) defined by 



Pi 



= XiP+, pt = XaP"^, Xi X2 = 1, 



Pi_L = xiPj_ A;_L, p2± = X2P1. - k±, 



(2.4) 
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the momentum-space wave-function "^^^^ can be expressed as 



^^^^(P1,P2, Ai, A2) = Rx^x.ix, k±) 



ss^ 



(2.5) 



where k±) describes the momentum distribution of the constituents in the bound state, 
and Rxfx2 constructs a state of definite spin {S, Sz) out of hght-front hehcity (Ai, A2) eigen- 
states. Exphcitly, 

RfX{x,k^) = Y.{\i\nij{l - x,k^,m,)\s,){\2\nUx,-k±,m2)\s2){ls,ls2\SSz), (2.6) 

where \si) are the usual Pauh spinor, and TZm is the Melosh transformation operator: 

rrii + XiMo + ia ■ kj_ x n 



7^A/(x, k±,mij 



^{rrii + x,Mo)2 + k] 
with n = (0, 0, 1), a unit vector in the z-direction, and 



2 + /c^ 1711 + ^ 



+ 



Xl 



X2 



In practice it is more convenient to use the covariant form for -Rff^^ 



V 2 Mq 



where 



Mo = - (mi - m2)2, 

r = 75 (pseudoscalar, 5 = 0), 
^ ■ (Pi - P2) 



Mo + mi + m2 



(vector, S = 1), 



with 



£^(±r 



-|^£l(±i)-Px, 0, £-i(±i) 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



Note that the longitudinal polarization 4- vector £^(0) given above is not exactly the same 
as that of the vector meson [cf. Eq. (|2.471 )]. We normalize the meson state as 



(M(P', S', S',)\M{P, S, Sz)) = 2{2iTfP+6\P' - P)6s'sh., 



(2.12) 
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so that 

dx (Pk I 



2(27r)= 



x,ki_)Y = l. (2.13) 



In principle, the momentum distribution amphtude 0(x, fc^) can be obtained by solv- 
ing the light-front QCD bound state equation |l5|JT7||. However, before such first-principles 
solutions are available, we would have to be contented with phenomenological amplitudes. 
One example that has been often used in the literature for heavy mesons is the so-called 
Bauer-Stech-Wirbel (BSW) amplitude |18[, which for a meson of mass M is given by 
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(x, /i;±)bsw = U^Jxil - x) exp { 1 exp 



(x - Xo) 



(2.14) 



2cj2 

where A/" is a normalization constant, x is the longitudinal momentum fraction carried by 

2 2 

the light antiquark, Xq = — -qj^)^ and u; is a parameter of order Aqcd- 
An other example is the Gaussian-type wave function. 



dk-^ I k"^ \ 

0(a;, A;_L)Gauss=A/'W-^exp -— , (2.15) 

where M = A^ir /tu'^Y^^, and k^ of the internal momentum k = {k±, k^) is defined through 

ei - 62 + k^ 
x = , l—x = , (2.16) 

ei + 62 Ci + 62 



with 6j = \/mf + k'^. We then have 



xMn m"^ + k'j , , 

Mo = 61 + 62, k, = ^f— ^, (2.17) 



2xMo 
and 

dk^ 6162 



(2.18) 



dx x{l — x)Mo 

is the Jacobian of transformation from {x, k±) to k. This wave function has been also used 
in many other studies of hadronic transitions. In particular, with appropriate parameters, 
it describes satisfactorily the pion elastic form factor up to ~ 10 GeV^ [Q. A variant of 
the Gaussian-type wave function is 



(..M-W§exp(-g). (2,19) 



with Mo being given by ( p.8|) . This amphtude is equivalent to 0(a;, A;_L)Gauss when the con- 
stituent quark masses are equal but becomes different otherwise. Nevertheless, we will not 
pursue this wave function further because it does not have an appropriate heavy-quark-limit 
behavior (see Sec. 111). 

Obviously, the Isgur-Wise function for heavy meson transitions depends on the heavy 
meson wave function (j){x, kj_) chosen. It turns out that in contrast to the Gaussian-type 
wave function, the BSW wave function fails to give a correct normalization for the Isgur- 
Wise function at zero recoil in P ^ V transition. 

A. Decay Constants 

The decay constant of a pseudoscalar meson P{qiq2) defined by (0|A^|P) = ifpP^ can 
be evaluated using the light-front wave function given by (2.1) and (2.5) 

(0|g27+75gi|P) = / {d?Pi}{d?P2}2{2Tif5{P - pi - P2)M^. k^)KxA^^ k^) 

X (0|g27+75gi|gig2). (2.20) 



Since Moyx(l — x) = y A"^ + kj_, it is straightforward to show that 

, -\/3 fdx(Pk± , , , , ^ 
fp= ——A 0^ fcx . 2.21 

V2J 2(27r)3 U2 + p 



where 

A = mix + m2{l - x). (2.22) 

Note that the factor y/S in ( p.21j ) arises from the color factor implicit in the meson wave 
function. 

Likewise, for the vector-meson decay constant defined by 

{0\V^\V) = fvMve^, (2.23) 

is found to be 

VS rdxd^k± (f)v{x,k^) I f ,^ . , ,2 , , ,2 

Tv = / — ; — , < xil — X Mt/ + rriimo + k , 

V2J 2(27r)3 .r^^^MvX ^ ^ ^ 1 2 X 



B 

+ 



2mi 



] ^ - 1 - 2x)Ml 

\ — X X 



(2.24) 



7 



where 



B = xnii - {1 - x)m2, Wy = Mq + rrii + 1712. (2.25) 

When the decay constant is known, it can be used to constrain the parameters of the hght- 
front wave function. 

B. Form Factors for P ^ P Transition 

With the hght-front wave functions given above, we will first calculate the form factors 
for P ^ P transitions given by 

(Psi^iPi) = uiq'm + + f-iq'm - ^2)^ (2.26) 

where = q2l^(li- For later purposes, it is also convenient to parametrize this matrix 
element in different forms: 



(P2\V\Pi) = v'A/iMjft+(,')(i.i + 1,2)" + - 1,2)" 



p. , p. ^ Ml_MiA ^ F,(f)Mt^,., (2.27) 



with Vi = Pi/ Mi, and 

Fi{q') = U{q% Fo{q') = U{q') + J^f^f-(^")- (2-28) 

In the heavy-quark limit Mi 2 — > 00, heavy-quark symmetry requires that [H 

h+{<l')=avi-V2), h^{q') = 0, (2.29) 

where ^(fi ■V2) is the universal Isgur-Wise function normalized to unity at the point of equal 
velocities: .^(1) = 1. The form factors Fi and Fq are related to the transition amplitude with 
the exchange of a vector (1~) and a scalar (0"*") boson in the t-channel, respectively. 

As explained in the Introduction, we shall work in the frame where g_L = so that 
= q~^q~ will cover the whole time-like region > 0. Define r = P2 / Pi (it is denoted by 
R in [|12I), then 

q^ = {l-r)(Ml-^\. (2.30) 



Consequently, for a given g^, there are two solutions for r: 



T± = ^[Vi-V2± J{vi ■ V2Y - 1 ) , (2.31) 



Ml 

where vi ■ f 2 is related to by 

The +(— ) signs in ( p.31|) correspond to the daughter meson recoiling in the positive (nega- 
tive) 2;-direction relative to the parent meson (call them the "+" and "— " reference frame, 
respectively). At zero recoil (g^ = (Jmax) maximum recoil (g^ = 0), r± are given by 
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r+(0) = l, r_(0) = (^) . (2.33) 

The form factors f±{q^) of course should be independent of the reference frame chosen for 
the moving direction of the daughter meson. For a given g^, suppose we obtain 

{P2\V+\Pi) =2P+i7(r+), (Psl^+I^i) =2P+i7(r_). (2.34) 

r=r+ r=r_ 

It follows from (|2.27| ) that 



2, (l-r_)iJ(r+)-(l-r+)i7(r_) 



r+ — r_ 

2, (l+r_)/7(r+)-(l + r+)i7(r_ 



/_(g^) = ' ^ . (2.35) 

r+ — r_ 

It is easily seen that f±{(f) are independent of the choice of "+" or "— " frame, as it should 
be. 

As noted earlier, in a frame with g+ > 0, there are actually two distinct contributions to 
the hadronic matrix element |^,T^,T^,T^]: valence (partonic) contribution calculated with rel- 



ativistic light-front bound-state wave functions, and non-valence (non-partonic) contribution 
(or the so-called Z-graph) arising from quark- ant iquark pair creation from the vacuum. In 
the following, we shall first provide some details for calculating the valence contribution, and 
then come back to the non- valence subprocess in Sec. II. D. For Pi = (gig) and P2 = (g2^), 
the relevant quark momentum variables are 
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p+ = (1 - s)P+, = xP^, p^^ = (1 - x)Pn_ + k±, = xPi± - fcj 

pt = {l-x')P^, p'^+ = x'P+, p,, = {l-x')P2± + k'^, p'^_^=x'P2±-k'^, (2.36) 



where x {x') is the momentum fraction carried by the spectator antiquark q in the initial 
(final) state. The spectator model requires that 



(2.37) 



Taking a Lorentz frame where Pi± = P2± = amounts to having q± = and k'j_ = k^. Then 
we readily obtain 



{P2\V^\Pi)= E /Kpj02(a;',M0i(a^,fcj 



Ai,A2,A 



00 



Substituting the covariant form given in Eq. into Eq. ( ^.381 ) yields 



{P2\V^\Pl) 



1 r dx (Pk I 



rJ 2(27r)3 



2Mo2MoiJ(l-a;')(l-a;) 



xTr 



75(^2 + "^2)7^(^1 + "^1)75(^5 - mq) 



After some manipulation, the trace term in the above expression is reduced to 



Tr 



where 



75(^2 + "^2)7^(^1 + "^1)75(^5 - rriq)] = -{A1A2 + kl)P+, 



Ai = rriix + rriqil — x), A2 = m2x' + mg{l — x'), 



and use of ( p.36|) has been made. Since 



MoiJx{l - x) MUx'{1 - X') = JAl + klJAl + kl 



A1A2 + kl 



we find from ( ^1341) , (^1381) , (|^ and O) that 

H(r) = dx — r:; , kj_)(j)i(x, k±) - 

^ Wo J 2(2vr)3'^2^ ^^'"'^ ^ ^Aj + kl^Al + kl' 
with x' = x/r. The form factors f±{q^) can then be obtained from ( |2.35|) . 



(2.38) 



(2.39) 



(2.40) 



(2.41) 



(2.42) 



(2.43) 
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As stated before, in the literature these form factors are customarily evaluated in the 
frame where = — P-^ = 0. This leads to = —qj_ < 0, implying a space-like mo- 
mentum transfer. The advantage of the condition g+ = is that form factors only receive 
valence contributions (see Sec. II. D). However, there are two drawbacks in this approach: 
First, form factors in the physical time-like region cannot be obtained without making ad- 
ditional extrapolation assumptions. Second, no information can be obtained for the form 
factor f~{q^) since P2 = Pi [see (2.26)]. At the maximum recoil q^ = 0, the form factor 
/+(0) is evaluated to be 

rf^/ci , AiA2 + kl 



(2.44) 



'Al + kl^Al + kl 

with Ai and A2 given by ( p.41| ) except for that x' = x here. Therefore, the results of ( p.35| ) 
and ( |2.43| ) at = are in agreement with ( |2.44| ). 



C. Form Factors for P ^ V Transition 



Form factors for P ^ V transition are defined as 

2 



{V{Pv,e)\J,\P{Pi)) 



Mp + M^ 



e-Pi 



ie,,^pe''P^PfV{q') - [(Mp + My)e,Ai{q 
e-Pi 



Mp + My 

where J^ = V^- A^, ^3(0) = Ao(0), 

2 Mp + My 



(Pi + Pv),A,{q^) - 2Mv-^ q, {A;{q') - A,{q^) 



and 



£'^(±11 



A,{q' 



2 el ■ Py, 
p+ 



2M 



V 



A.{q' 



Mp - Mv 



2M 



V 



Mq' 



-Ml + P 



VI- JD+ JD 



! -< y ; 



(2.45) 



(2.46) 



(2.47) 



are, respectively, the transverse and longitudinal polarization vectors of the vector meson. 
The form factors Ai and A2 are related to 1"*" intermediate states, Aq to 0^ states, and V to 
1^ states. The P —* V matrix element can also be parametrized in different ways: 

{V\J^\P) = ig{q^)e,,^pe''Pi^P^ + f{q^)e^ + {e ■ P^)[a^{q^){P^ + Py)^ + a_{q^){P, - Py)^] 

lMpMv{ z^(g') W^'^'"^^ + h<fy^^ 



+ {e ■ v)[a+{q^){v + v')^ + a_(g')(t; - v')^] 



(2.48) 
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where v = Pi/Mp and v' = Pv/My- They are useful for later discussions. The form factors 
a±, f and g are related to V, v4o,i,2,3 via 

2 



Mp + Mv 
1 

Mp + My 

2Mv 



-{Mp + Mv)Ai{q' 



q' 

2Mi 



As{q') - Ao{q')] 



V 



Mp + My 

2Mv 



A,iq' 



Mp - Mv 



2M 



V 



A2{q')-Aoiq' 



Also, 



Mp + Mvg{q^). f{q') 



fiq' 



^Mp + M 



V 



Ml 



(a+ + a_), 



(2.49) 



(2.50) 



MpMv (a+ - a_). 



VMpMy 

In the heavy-quark limit Mp, My — oo, heavy-quark symmetry demands that 



a+ + a_ = 0, 



~g = av-v'), f = -{l+vv')av-v'). 



(2.51) 



The calculation of the P ^ V form factors is more subtle than the P ^ P case. If we 
choose a frame where Pi± = Pv± = as before, we will have e ■ Pi = for transverse polar- 



ization. As a result, form factors a± in ( |2.48| ) cannot be separately determined. Therefore, 
we will let Pi± = Pv±^0 at the outset, and set them to zero only after the form factors are 
extracted. With the transverse polarization £^(±1), form factors a±{q'^) and giq"^) can be 
individually determined. Then using the longitudinal polarization e^{0), we are able to fix 
the remaining form factor fiq"^)- 

We begin with a±(g2). Since £+(±1) = [cf. (gTzl) ], it follows from (|]48D that 



- {V{Py)\A+\P{P,)) = [e ■ P,)[a^{Pt + P^) + a_{Pt - P^)] 

(^"1 ■ P±)[a+(1 + r) + a_(l - r)]P+, (2.52) 

with r = Py /Pi and Pi = Pn = Pyj,. As will be shown below, the above matrix element 
at the quark level has the form 



{V\A-^\P) = 2e± ■ Pi(l - r)/(r)P+ 



(2.53) 
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Substituting this into (|2.52|) and solving the equations for r = r+ and r = r_ yields 

2s ^^+(1 - r-)I{r+) - r_(l - r+)/(r_) 
a+{q ) - 



r+(l + r_)/(r_,_) — r_(l + r_,_)/(r_) 



(2.54) 

r+ — r_ 

in analog to Eq.( |2.35| ) for f±{q^). In order to illustrate several subtle points in the deriva- 
tion of /(r), we will go through the calculation in a bit more details. First of all, it is 
straightforward to show that for P = {qiq) and V = {q2q) 

dx<fk\ —2x' 



{V\A+\P) 



2(27r) 



lAl + kl^Al + k'l 
where Ap = Ai, Ay = A2 [see Eq. (|2.41| )], and 



(f)y{x , k^)(f)p{x, k±){a + &), 



(2.55) 



a = mi{e- + i ■ PqVt) + ^2{s ■ PiPq - e ■ p^^pf) + mg{i ■ p2pt + i ■ Pipt), 

, e - {p2-Pg) , + + + + + + 

b = — {minigpj - m2mqPl - mimapj + Pi ■ PgPj - Pi ■ P2PJ + P2 ■ PgPi 

Wv 



(2.56) 



with e'^ = £'^(±1) given by ( p.47| ) and Wy = Mqv + 1712 + rriq. By virtue of (|2.36|) we find 
that 



a 
b 



[l-r){l-2x')Ap{e^-P^)P, 
ApBv + k^ 



-2 1 - r 



Wv 



s± ■ P±)P^ + 



with 



By = —m2x' + (1 — x)mq. 
The ellipses in ( ^.571 ) denote contributions from terms proportional to el ■ /c^ in 



(2.57) 



(2.58) 



Naively, these terms linear in are not expected to make contributions after integrating 
over k 1 . But this is not the case. Consider the term 



(t)v{x\k'^ 



Al + 



(2.59) 



and note that fc^ is different from k± due to a non- vanishing Pj^: 

k'^ = + (x' — a;)Pi, 



(2.60) 



where we have used (|2.36|) and ( |2.37|) . Consequently, 
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V(fc^) = Mkl) + {d(l)y/dkl){k'l -kl) + -- 



(2.61) 



with 



Since 



(pv 



dki 



(2.62) 



d'k^ie^ ■ k^){k^ . P^) = ^ / dPk^kl{e^ ■ Pj_ 



(2.63) 



it is evident that the hnear term {e±_ ■ k±) in ( |2.57|) will combine with the linear term {k± ■ P_ 



in ( p.61|) to make a contribution to (V^jA+lP). We wish to stress that this additional contri- 
bution from 0v was first noticed and obtained by O'Donnell and Xu [0,0 and was neglected 
in the work of Jaus [|6|,p!0|. 

By the same token, in the expression of 6, the {e±_ ■ k±) term in i ■ {p2 — Pq) will also 
combine with the {k± ■ P±) term in {mimqP2 + ■ ■ ■) to yield a contribution proportional to 
e± ■ P± after integration over k±. The final result is 



/(r) 



dx I ^^'^^^^^'^^^'^^^^'^^^ I (1 - 2x')Ap 



{ApBv + kl){l + ev{x',k^)kl) + \e 
Wv 



+ [(1 - 2x')Ap - Av]Qv{x\ kA_)kl - 2 



^2.64) 



with x' = x/r. In deriving (|2.64| ) we have first integrated out x' and k'^. We can of course 
alter the order of integration by first integrating over x and k± and obtain 



/-i /• d^k^ x(j)*vix',k±)(j)p{x,k 



Al + klJAl + kl 



X \Av-[{1- 2x')Ap - Av]&p{x, ki_)k\ + 2 



{ApBv + kl)Qp{x, k^)kl + \k^ 



V 



(2.65) 



with X = x'r, where we have used the notation k± instead of k'j^ as it is a dummy variable. 
The result ( |2.65| ) will be utilized in Sec. 3 to show that light-front model calculations fulfill 
the heavy-quark-symmetry requirement ( |2.51| ). 

At g2 = 0, we have r+ = 1 and r_ = {My/Mpf. It follows from (|23|) and (g]^ that 
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^2(0) = (Mp + My)a+(0) = -{Mp + Mv)I{r = 1) 



"^^y 2(27r^3 



d'^k± X(j)Y{x, k±)(j)p{x, k±) 



(Mp + Mi/)Ul-2x)^P 



(2.66) 



+ [(1 - 2x)^p - ^y]0v(x, A;x)A;2 - 2 



V 



This is in agreement with Eq.(30) of 0, but disagrees with the result obtained by Jaus P|. 

Having fixed a±(g^), we are ready to calculate f{q^) in ( p.48|) . From (|2.49| ) it is clear 
that once f{q^) is determined, so are the form factors Ai{q'^) and AQ{q^). Since the "+" 
component of is needed to extract f{q^)i we consider the longitudinal polarization £^(0) 
of the vector meson V and take a frame where Pi = 0. Hence, 



-{V\A-^\P)= f{q')—P^ 

Let {V\A+\P) = J{q^)P^, then 

1 



2M,. 



V 



[a+(l + r) + a_(l-r)]P+ 



Mi 



[a+(l + r) + a_(l-r)]-^J(g2). 

r 



2 \ - r ^ 
After a straightforward manipulation, we obtain 

d'^k± x'(f)y{x',k±)(f)p{x,k 



-r / dx 
Jo 



2(27r) 



(c + rf). 



with 



1 



(l-x')-MoV^P + 



1 



Mov xWv 
and [cf. Eq.(pl)] 



(ApBv + k] 



Ap + k,(mi H mq) 

X X 



ml + kl ml + kj 



X' 



X 



ml + k\ ^ml + k]_ 



1 — x' ' x' 

To check the above results, we note that for r+(0) = 1 [see (|2.33|) ] 



(2.67) 



(2.68) 



(2.69) 



(2.70) 



(2.71) 



{V\A'^\P) 



r(0)=r+(0) 



-[/(O) + (Ml - M2)a+(0)]P+/Mv = 2Ao(0)P+, 



so that 



(2.72) 



^(0) = 2"^(0) 



r=l 



(2.73) 
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where use has been made of ( p.49|) . Then it is not difficult to show from ( p.69[ - pl70| ) that 

dxd^k^ ApAy + (1 - 2x)kl + H(!^i±!^ 



^(0) 



Wv 



2(27r)3 



Al + fc2 ,IAI + fc2 



(2.74) 



which agrees with Eq.(27) of ^ obtained in the = frame (implying r = 1). It will 
be shown in Sec. III.B that our result for f{q^) does respect the heavy-quark-symmetry 
requirement. 

Thus far we have imposed the condition = to extract the form factors a±{q^) and 
f{q^)- For the vector form factor g{q'^) or V(g^), it proves more convenient to first let g_L 7^ 
and then set it to zero after the vector form factor is obtained. The "+" component of the 
vector matrix element for transverse polarization reads 



= zge+^,y[-e-q^py + (Pi - PyYe^py - P{e^qy 



\x 



(2.75) 



where Pj_ = Pi_l, Pyx 
{V\V+\P) 



P± — q±. At the quark level, we have 

d^k± 2x'4)*y{x', k±)(f)p{x, ki_) \e ■ {p2 - Vq) 



J 2{2txY 



Al + klJAl + kl 



V 



(2.76) 



The transverse momentum variables are 



Pii. 



'l-x)P^ + k, 



P2± 



[1 - x)P± - g_L + k^, pg± = xP_L - k_i 



(2.77) 



It suffices to set = e^{±) in ( p.76| ) to get contributions proportional to e+-xye~{q^py — 
q^p^), which is related to the ffist term of ( |2.75|) at the hadron level. It is easy to check that 
the transverse components e^(±) and £^(±) will not generate the same structure. Repeating 
the similar derivation as before, we obtain 

2V{q^) r, f d^k^ 2x'(p*y{x',k±)^p{x,k^) 



9[Q 



Mp + Mv 



"^""J 2(27r)3 



Al + klJAl + kl 



X \Ap + {Ap- Ay)ey{x', k^)kl + 



rkl + (1 - r) l^xMopk^ 



'U2 



X k 



l-x' 



(2.78) 



+ (1 - r)Qy{x', ki_)ki{2x^Mlp - x'^M^y - ml - kl) 
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with being defined in ( CT )- For r(0) = r+(0) = 1, (|2T8|) leads to 

X l^^p + ^ + - m2)ey(x, fc±)A;ij , (2.79) 

which agrees with 

Therefore, we have calculated the form factors f{q^)i g{q'^) and a±(g^) in the time-like 
region within the light-front framework. Form factors Vlq"^) and Aq, 1,2(0'^) can then be 
determined via Eq. (|2.49|) .P] 

D. Non- valence Contribution 

Thus far we have concentrated on the valence-quark contribution to the form factors. As 
stated in the Introduction, there also exist contributions which are generated from the quark- 
antiquark excitation or higher Fock-states in the hadronic bound states. This additional Z- 
graph contribution vanishes in the frame where the momentum transfer is purely transverse 
i.e., = 0, but survives otherwise. 

The general feature of the non-valence configuration can be recognized by considering 
the quark triangle diagram (see Fig. 1). In terms of the "+" component of momenta, 
the Feynman triangle diagram in the light-front framework consists of two subprocesses: 
one corresponds to the valence-quark approximation for the meson wave functions, and the 
other to the contribution of quark-pair creation from the vacuum. That is, through the 
mechanism of quark-antiquark pair creation, the "spectator" quark in the second subprocess 
is fragmented into a meson plus an outgoing quark. A detailed study of the quark triangle 
diagram for P ^ P transition gives (generalization to P ^ V transition is straightforward) 



{P2\q2l^qi\Pi) = Mt + Mt, (2.80) 



^In the frame where = 0, only three of the P ^ V form factors, namely /, g and a+ or V, Ai 
and A2 are determined. However, Aq can be fixed at = using the relation ^o(O) = ^3(0) and 

dH). 
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with 



= '''' lo x(l - x)(l - X') I 



dx f d'^ki N;! 



lo x{l-x){l-x')J 2(27r)3 {M? - Ml,){Ml - Mi,y 

^+ ^ _ dx r d^k± iVj 

^'^'Jr x{l - x){l - X') J 2(27r)3 {Ml - Mi,){q^ - Mh)^; ^^"^'^ 

where x' = x / r , gi and g2 are the quark-meson couphng constants at different vertices, Mi 
and M2 are the masses of the initial and final meson respectively, Mq^ (^02) is the same as 
M^p {Mly) defined in (gTTD , and 



\ \ — X r — X J 

= A[pl{m2niq + p2 ■ Pq) + ptimiTrig + pi ■ pg) + p^{mim2 - pi ■ P2)], (2.82) 

= 4[p^{-r]Mim2 + Pi ■ P2) +p^(-r7Mimi + Pi ■ pi) + P^{mim2 - Pi ■ P2)], 

with 7] = {mi — mq)/Mi. Since M.a receives contributions from the kinematic region < 
X < r or < /c^ < P2 (see Fig. 1), it corresponds to the valence-quark configuration. As 
for A^^, only the region r < x < 1 or P2 < < Pi is relevant, and it corresponds to the 
non- valence contribution. It is straightforward to check that, apart from a sign difference, 
A''^ is precisely the trace term given in ( [2.401 ). This implies that the previous calculation 



for P —* P form factors in the Hamiltonian light-front approach is identical to the Feynman 
triangle graph under the valence-quark approximation. Obviously, making the following 
substitutions 

V^gi 1 4>i{x,k_L) 



x{l - x) Ml - Ml p^^^ 
V2g2 1 Mx',k^) 



x'{l - X') Mi - MI2 .L^+k 



(2.83) 



in M+ will reproduce the resuh {P2\q2l^qi\Pi) = 2P^H [see (pTS^ ) and ( 121431) ]. 

Unlike the valence-quark contribution, only the wave function (f)i{x,k±) of the initial 
meson enters into the expression of A^^; (f)2{x', k±) is not applicable for the non-valence graph 
because the light-front momentum of the spectator quark is larger than the momentum 
P^ of the daughter meson (see Fig. 1). This makes the task of calculating the effect of 
the Z'-graph considerably more difficult. Nevertheless, some qualitative features of Ai'^ 
can still be comprehended. First of all, as noted earlier, the contribution from non-valence 



configurations vanishes in a frame wliere g"*" = or r = 1. However, tliis frame is suitable 
only for space-like g^. Second, it is easy to show that — in the limit of heavy quark 
symmetry mq oo, because it takes an infinite amount of energy to create a heavy quark- 
antiquark pair. This has the important implication that we do not have to worry about the 
pair-creation subprocess when calculating the Isgur-Wise function. Beyond the heavy-quark 
limit, it is commonly argued that the non- valence contribution leads to a small correction in 
heavy-to-heavy transition but becomes more important for heavy-to-light decays |]l5| , |T0| , pr3 



For example, a i?*-pole contribution is usually believed to be the dominant non-valence effect 
in 5 — s> TT transition, especially when is near the zero- recoil point ||2^. Some estimates 
based on the i?*-pole contribution with the help of chiral perturbation theory indicate that 
for large values of g^, the Z-graph provides the dominant contribution to i? — > vr form factors 
(for a recent estimate, see [T^). 

In this paper we will demonstrate that even for heavy-to-heavy transition, the importance 
of the non- valence contribution depends on the recoiling direction of the daughter meson. As 
shown in (|2.31|) , for a given g^ there are two possible reference frames characterized by ^^(g^) 
and r_(g^), corresponding to whether one chooses the velocity of the final meson to be in 
the positive or negative z-direction relative to the initial meson. Of course, the form factors 
are independent of the choice of the "+" or "— " frame. This means that the combination 
M.'^ ( |2.8CI| ) should be independent of the choice of r(g^) = r+(g^) or r(g^) = r_{q^). 

From ( |2.81| ) we see that the r dependence of or A^^J" appears both in the integrand and 
in the integration limit. As a consequence, M.'^ and M.^ separately are in general "±" -frame 
dependent. In other words, the valence-quark and non-valence contributions to mesonic form 
factors are in general dependent on the recoiling direction of the final meson, hut their sum 
is not. For the form factors f±{q^) in P —>■ P decay and a±(g^) in P ^ \^ transition, we have 
"demanded" that the valence contribution itself be frame independent [see ( p.34|) , ( p.35| ) and 
( p.54| )]. For form factors Aq, Ai and V in P —>■ V decay, explicit calculations in Sec. IV. B 
show that the valence contributions for r = r_|_ and r = r_ are indeed different (see Fig. 6). 
Thus in principle we cannot make firm predictions for these form factors even for B ^ D* 
transition, unless the non- valence contributions are also calculated. Nevertheless, corrections 
due to the non-valence configuration are expected to be marginal for heavy-to-heavy form 
factors evaluated in the frame where r = r_|_, but become more significant in the "— " 
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frame (r = r_). The argument goes as follows: we know that the non- valence contribution 
vanishes if = 0. Now g+ is never zero in the "— " frame, whereas in the "+" frame g"*" = 
when r+ = 1 [see ( p.33| )]. That means the valence-quark contribution in the "+" frame is 
exact at the = point. As will be shown in Sec. IV. B, the valence contributions at = 
in the "— " frame are generally smaller than those in the "+" frame; the difference should 
be accounted for by the non-valence configuration. These points will be elucidated in more 
detail in Sec. IV.B. 



III. THE ISGUR-WISE FUNCTION 

In Sec. 2 we have computed the P P form factors f±{q'^) and P ^ V form factors 
V{q^), y4o,i,2(g^) for the entire physical q^ region using the light-front wave functions. It is 
very important to check if the light-front model predictions are in accord with the require- 
ments of heavy-quark symmetry, namely ( 2.29 ) and ( p.51 ). In other words, as niq oo, we 



would like to see if there exists a universal Isgur-Wise function which governs all heavy-to- 
heavy mesonic from factors in the light-front quark model. 

To our knowledge, the Isgur-Wise function has not been calculated directly for q^ > 
within the framework of the light-front quark model, though it has been considered in 
||9| JT9|JTT|| . The analysis of ^ is based on the observation that the knowledge oi P ^ P 



oi P -^V form factors at = (or at any point of q^) suffices to determine the Isgur-Wise 
function in the whole kinematic region. However, this relies on the assumption that the 
model calculations of form factors obey heavy-quark symmetry and that the universal form 
factor is only a function of v ■ v'. The Isgur-Wise function is derived in |jT^ from space-like 
elastic form factors of heavy mesons,^ while it is obtained in [jlT| by performing an analytic 



continuation from the region < to time-like momentum transfers. In contrast, we do not 
impose heavy-quark symmetry from the outset, so that we can check explicitly if the weak 
decay form factors of heavy mesons can indeed be described by a single universal function 
when rnq oo. We will calculate this universal function directly at the time-like momentum 



^This is based on the argument that, for the elastic form factor, q'^ = —{y-v' — 1)/(2M). Thus the 
space-Uke elastic form factor is related to the Isgur-Wise function at time-like momentum transfers 
{vv' > 1). 
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transfer to see if it is independent of heavy quark masses and their ratio. It is important to 
note that, since heavy quark-pair creation is forbidden in the mq —>■ oo hmit, the Z-graph is 
no longer a problem in the reference frame where g+ > 0. Therefore, within the light-front 
quark model, we are able to compute the Isgur-Wise function exactly for time-like g^. 

To proceed, we first investigate the heavy- quark- limit behavior of the wave function. In 
the infinite quark mass limit mq oo, the light-front wave function has the scaling behavior 
[17]: 

(pQg-{x,k^)^ ,/m^^{mQX,k^), (3.1) 



where the factor ^ttiq or vM (M being the mass of the heavy meson) comes from the 
particular normalization we have assumed for the physical state in ( p.l2[ - pTT3D . The reason 



why the light-front heavy-meson wave function should have such an asymptotic form is as 
follows. Since x is the longitudinal momentum fraction carried by the light antiquark, the 
meson wave function should be sharply peaked near x ~ Aqcd/'"^q- It is thus clear that 
only terms of the form "rngx" survive in the wave function as mg — oo; that is, tuqx is 
independent of tuq in the mg — > oo limit. For the BSW wave function (|2.14|) , we find that 

$(X, /.x)bsw = f ^) exp exp i-^ (3.2) 



where X = rnqx, and the normalization condition ( p.l3| ) becomes 



/o J 2(27r)3 
For the Gaussian- type wave function ( p.l5| ), we obtain 



iX —^AHX,ki_)i' = l. (3.3) 




*(X,*Jo.„„=4(-^)''\xp -il)^=|. (3.4) 



From ( |2.17| ) it is clear that = [X — (m? -|- kY)/X]/2 in the heavy-quark limit. Therefore, 




However, since Mq — * itlq + 0{mQx), it is clear that the wave function (|2.19|) , which is a 
variant of the Gaussian type, does not have the correct asymptotic form in the heavy-quark 
hmit. Hence it is not suitable for describing heavy-quark transitions. 
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A. P ^ P Transition in Heavy-Quark Limit 



With the hght-front wave function k^) constructed in the rriQ —>■ oo hmit, the P — 
P transition function H{r) in the hmit of heavy-quark symmetry (i.e., mi,m2 oo 

becomes 

AiX)AiX') + kl 



<Pk, 



H(r) = J— / dX , , 
^ ' ^ MiJo J 2(27r)3 



$(X',fc^)$(X, k^) 



A''{X) + klJA''{X') + k 



(3.6) 



where X = rriix, X' = m2x' , and A{X) = X + mq. Note that the quantities X, X', rriq 
and A;_L appearing in the integrand are aU of order Aqcd- Denote z = V2 /vi = {Mi/M2)r, 



then [see (|23l| )] 



Z±=Vi-V2± \/iVi ■ - 1 . 



(3.7) 



Obviously, z+z^ = 1 and X'/X = l/z. Let H{r±) = y'Ma/Mi H{z±), so that (|]3|) can be 
rewritten as 



{P2\V+\P,) = 2JMiM2H{z)vt. 



By a simple change of integration variable, one can readily show that 



(3.^ 



H{z) = zH{l/z). 



(3.9) 



To check the validity of the heavy-quark-symmetry relation ( p.29| ), we note that h±{ 



are related to H{z) via 



l-z^ 



H{z)±zH{l/z) 



(3.10) 



in analog to ( |2.35| ) for f±{q^). By virtue of (PT^), the HQS relation /i-(g^) = given in ( f^.29D 
is indeed satisfied, and the Isgur-Wise function is given by 

2H{z) 



i{vi ■ V2) 



l + z 



(3.11) 



Evidently, the Isgur-Wise function is independent of the heavy quark masses mi, m2 and their 
ratio, but it depends on the light spectator quark mass. The R.H.S. of (|3.11 ) is invariant 
under the exchange z ^ l/z, implying that the Isgur-Wise function ^(t>i ■V2) is independent 
of the choice of the recoiling direction of the daughter meson, as it should be. At zero recoil 
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{z = 1), the expression for H{1) becomes identical to the normahzation condition (|3.3|) . 
Hence H{1) = 1, and the Isgur-Wise function obeys the correct normahzation condition 
^(1) = 1. We would like to stress again that, unlike the previous works 0,|lT| where ^{v ■ v') 



is actually evaluated for B —>■ D transition and for space-like values of g^, here the Isgur- 
Wise function is obtained in the infinite quark mass limit and calculated directly for > 0. 
Within the specific model we have taken, our result is exact. 

In the limit of heavy-quark symmetry, form factors Fi and Fq are related to the Isgur- 
Wise function via 

^(^1 ■ ^2) = , Fi{q)= TT-TT T. 7^2 T • (3-12) 



M1 + M2 Ml + M2 



1 - 



(Mi+M2)2_ 

Hence the dependence of Fi is different from that of Fq by an additional pole factor. 

B. P ^ y Transition in Heavy-Quark Limit 



There are four HQS relations given in ( |2.51D for P ^ V form factors. We shall first focus 



on a± (or a±). As 1711,1712 — > 00, we can show that 

Qvix', k^) ^ e(X', k^), ep(x, k^) e(X, A;^), (3.13) 

with Qy being defined in ( p.61| ). The P V transition amplitude /(r) [see ( p.64| ) and 
( p.65| )] reduce to 

I{r) = I{z) = -^^^ J ^(^[^'^(^) + ^'(^ - ^')0(^', k^)kl]f{X)f{X') 
dX'cPk 



^mim2 J 2(2 



/ ^fl^f^^^^'^ " ~ (3.14) 



where /(X) = $(X, ki_)/ ^ A^{X) + k\, X' = X/z, and all terms proportional to l/Wy have 
been neglected in the heavy-quark limit. It is evident that I{z) satisfies the relation 

I{z) = I{l/z). (3.15) 

Therefore, from ( |2.54| ), 

2Ml r+r_I(z+) - r+rJ{z\ 
a+ + a-= .^^ — ^-^ ^ = 0, 3.16 
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and 

a+ - a_ = 2^MpMv '"'^'^^^'^^ ^^'^^^'^ = ~2^MpMvI{z). (3.17) 
By comparing this with ( ^.511 ) yields the Isgur-Wise function 

C(v.v') = 2rdX ^{X',k,)^XM) 

^0 J 2(27r)3 ^A''{X) + kl^A\X') + kl 

X [X'A{X) + X\X - X')Q{X\ k^)k]_] . (3.18) 

with X' jX = 1/z. It remains to show that C,{v ■ v') is indeed the same as the Isgur-Wise 
function ^(t> ■ v') found in P — P transition ( p.ll| ). We will address this issue later in Sec. 
IV. After showing the HQS relations ( p.l6| ) and ( p.l7| ) for form factors a±, we turn to the 



vector form factor. One can easily show from (|2.78| ) that, indeed, 



~g{q^) = ^MpMygiq') ""^'^ -2^MpMyI{z) = ({vv'), (3.19) 
in accord with ( |2.51|) . 

Using the results ( p.l6| ) and ( p.l7| ) for form factors a±, we are ready to prove the remaining 



HQS relation for f{q^). It follows from (^ABj ), (|23|) and (g:69| - gT0D that 



f( 2^ 1 X'~X'' r dXd'k^ 1>(X^A:x)$(X,fc^) f My \ 

■^^^^ v/MWl^ 2 XX' ^+7 2(2vr)3 ^J^^^W^^J^^^K ' ^ " ^ 

where terms proportional to 1/Wy vanish in the limit of heavy-quark symmetry. We find 
from (|270|) that 

x'My HQ limit , , ,2n /o r,-, \ 

-c ^ -2—{ApAy + kl), (3.21) 

TTll A 

hence 

1 v2 Y' / X \ 

he) = "2 ^3^c(- • - X + x^j ■ ^^-^^^ 

where use of (|3.11| ) has been made. Then, using X'/X = 1/z and (^3)> 



desired HQS relation given in ( p.51|) : 



f{q^) = -(l + vv')av-v'), (3.23) 
provided that ({v ■ v') = i{v ■ v'). 
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Is the function ({v ■ v') given by (|3.18|) identical to the Isgur-Wise function ^{v ■ v') ? 



? 



While ^(1) = 1 is always valid irrespective of the details of the light-front amplitude used, 
the normalization of ({v ■ v') at zero recoil is nontrivial. In fact, we find that ({1) depends 
on the choice of the light-front model wave function. We find numerically (see Sec. IV.A) 
that the HQS requirement ({1) = 1 is fulfilled by the Gaussian- type wave function (|2.15| ), 



but not so by the BSW-type wave function ( p. 14 ). In other words, the normalization of 



the Isgur-Wise function at zero recoil in P ^ V transition puts a severe restriction on the 
phenomenological light-front wave functions. Since we are not able to solve the light-front 
QCD bound-state equation to obtain the momentum distribution amplitude k±), we see 
that heavy-quark symmetry is helpful in discriminating between different phenomenological 
amplitudes. As will be shown in Sec. IV.A, ({v ■ v') is numerically equal to ^(f ■ v') if the 
Gaussian-type wave function is used. 

The P ^ V form factors in the heavy-quark limit are all related to the Isgur-Wise 
function via 



- Mp + My ^^^^ ) - Mp + Mv ~ ^ • ^^-^^^ 



1 - 



That means V, Aq, A2 all have the same dependence and they differ from Ai by an 
additional pole factor. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

To examine numerically the form factors derived in the last section, we need to specify 
the parameters appearing in the phenomenological light-front wave functions. We shall use 
the decay constants to constrain the quark mass and the scale parameter u. The decay 
constants of light pseudoscalar and vector mesons areQ 

/^=132MeV, /K=160MeV, /p=216MeV, /x*=210MeV. (4.1) 



^The decay constant fp is obtained from the measured decay rate of e^e , while fx* is 

determined from r K*Ut-- 
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The decay constants of heavy mesons are unknown experimentally, so we have to rely on 
model calculations and lattice results. To be specific, we take 



fo = 200 MeV, /b = 185 MeV, fo* = 250 MeV, fe* = 205 MeV, 



(4.2) 



where the estimates for Jd* and Jb* are relatively more uncertain. The parameters rUg and u 
in the Gaussian-type and BSW-type wave functions fitted to the decay constants via ( p.21|) 



and ( 2.24 ) are listed in Table I. Note that the quark masses given in Table I are fixed to 
the commonly used values, and the other fitted values are by no means unique. Presumably, 
other hadronic properties, for example the light-meson elastic form factor measured at a 
wide range of momentum transfer, would be helpful in fixing the light-front parameters. 



Table I. Parameters niq (in units of GeV) and u! in the Gaussian-type and BSW-type wave 
functions fitted to the decay constants given by ([4.1D and ([4.2|). 



wave function 






rric ujd l^D' 




Gaussian 
BSW 


0.25 0.33 0.30 
0.25 0.30 0.29 


0.40 0.38 0.31 
0.40 0.34 0.30 


1.6 0.46 0.47 
1.6 0.46 0.46 


4.8 0.55 0.55 
4.8 0.58 0.57 



A. Results for the Isgur-Wise function 



Before proceeding to numerically evaluate the P P and P ^ V form factors, it is 
important to check the Isgur-Wise function to ensure that model calculations do respect 



heavy-quark symmetry in the infinite quark mass limit. With the Gaussian-type ( |3.4| ) and 
BSW-type ( p.2| ) wave functions given in the limit of heavy-quark symmetry, the Isgur-Wise 
function C,{v ■ v') for P — > P transition calculated from ( |3.11| ) and ( |3.10| ) is shown in Fig. 2 
using ojd = ^^b = 0.55. We see that the Isgur-Wise function obtained from Gaussian-type 
and BSW wave functions is very similar. The slope of ^{v ■ v') at the zero-recoil point is 



-e'(i) 



1.24. 



(4.3) 



Recent theoretical estimates and experimental analyses favor < 1. The slope parameter 
is subject to constraints from Bjorken and Voloshin sum rules (for a review, see 0). A tight 
bound is derived to be 0.5 < < 0.8 [21]. QCD sum-rule results range from 0.70 to 1.00 
21| . It thus appears that our slope parameter ( [4.3|) is too large. This may be attributed to 
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the fact that the Gaussian-type amphtude does not have enough amount of high-momentum 



components at large k±. It has been shown in that the one-gluon-exchange interaction 
can generate high-momentum components in the meson wave function and reduce the value 
of significantly. 

Although ^ is independent of heavy quark masses, it is interesting to see if it can be fitted 
to a simple pole behavior for a specific transition, e.g., B D : 

^^^^ (1 - gVMp^ole)" ' ^^-^^ 

where v ■ v' = (M| + M|, - q"^) /{2MbMd). We find that ^(g^) is fitted very well over the 
entire > region with a dipole behavior (one cannot tell the difference between fitted and 
calculated curves) with 

a = 2, Mpoie = 6.65 GeV. (4.5) 

Indeed, this pole mass is close to the mass 6.34 GeV of the l"*" vector meson with (be) content. 

The most interesting and striking results are shown in Fig. 3 for the function ({v ■ v') [see 
( ^.18| )] for P V transition obtained by taking the heavy-quark limit of the form factor 
^2(0'^) or V{q'^). For the Gaussian-type wave function, we find that ({1) = 1 at zero recoil, 
and that numerically ({v ■ v') is identical to ^(f ■ v')^ In contrast, the curve computed using 
the BSW amplitude deviates consistently from ^{v ■ v'); in particular, C(l) = 0.87 at zero- 
recoil. That means, for reasons not clear to us, the overlapping of the BSW wave functions 
for P ^ V transition at zero recoil is not complete in the heavy-quark limit. This in turn 
implies that the light-front amplitude $bsw is inconsistent with heavy-quark symmetry for 
P V transition. 

We note that the presence of the B term in ( |3.18| ) is crucial for obtaining the numerical 



equivalence of ({v ■ v') and ^{v ■ v'). Hence the form factors V^q"^), Ai{q'^), ^2(5^) obtained 
previously in [|6|,p!0| are incomplete since the G terms are not taken into account there. 



^Since numerically C{v • v') is equal to ^{v ■ v') up to six digits for the Gaussian-type amplitude, 
we believe that this equivalence is exact, although both Maple and Methematica fail to give an 
analytic result for (p.l8|) . 
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B. P ^ P Form Factors 

Since the BSW wave function fails to give a correct normalization at zero recoil for the 
Isgur-Wise function in P ^ V transition, the ensuing calculations are all carried out using 
the Gaussian-type wave function. The dependence of the form factors Fi{q^) = f+{q^) 
and -Fo(q'^) for B ^ D weak transition computed using ( p.28|) , (|2.35|) and ( |2.43|) are shown 



in Fig. 4 (we have neglected the non-valence contributions). At = 0, we obtain Ff^{0) = 
F(f^(0) = 0.70 . From Fig. 4 we see that F^^{q^) can be fitted by a dipole approximation 
in the entire time-like q^ region with a pole mass M^^Xe = 6.59 GeV, in agreement with the 
pole mass 6.65 GeV fitted to the Isgur-Wise function [cf. (|45| )], while F^^^q^) at low q^ 
(0 < ^ 6 GeV^) exhibits a monopole depenfence with Mpoie = 7.90 GeV. This monopole 
behavior for Fq^ at low q^ is consistent with (|3.12|) . 

The dependence of the form factor f+{q^) for the transitions B ^ n, B —>■ K, D —>■ it 
and D ^ K are shown in Fig. 5. The numerical results for the form factors at = are 

/^-(O) = 0.29 , /+^(0) = 0.34 , /:^"(0) = 0.64 , /+ ^(0) = 0.75 . (4.6) 

From Fig. 5 we see that, near the zero- recoil point, the valence-quark prediction for f^"^ 
decreases as increases. As explained in |T2|, the dipping of the valence-quark contribution 



toward the q^^-^ point can be understood as follows. Recall that the decay amplitude involves 
an overlapping integral of the wave functions of the initial and final mesons. If both mesons 
were heavy, then it is obvious that, by heavy-quark symmetry, maximum overlapping must 
occur at the zero-recoil point. However, in the situation of -B — * tt transition, the internal 
momentum distributions of the heavy B meson and light pion peak at different values of 
X. Specifically, 4'B{x,k±) has a narrow peak near a; = 0, whereas 07r(a;, fc^) peaks with a 
much larger width at x = 1/2. Consequently, maximum overlapping of the wave functions 
actually occurs somewhat away from the zero-recoil kinematics. For D it transition, 
maximum overlapping occurs in the close vicinity of zero recoil (see Fig. 5). Since the non- 
valence contribution is expected to be important for heavy-to-light form factors, especially for 
B IT transition, a comparison with data at large cannot be made until such contribution 
is included (form factors at g^ = are not affected by the pair-creation configuration) 



^In ||Tl| form factors at < are reformulated as a double dispersion integral representation, 
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We see from ( [4 .61) that while the predicted /+^(0) is in nice agreement with experiment, 
/+^(0)expt = 0.75±0.03 |2|; /^'^(O) and the ratio R = /^""{O)/ /^^(O) = 0.87 are too small 
compared to the measured values R = 1.29±0.21±0.11 and 1.01±0.20±0.07 |2|]. It has 
been pointed out in that the unexpected large decay rates of Cabibbo-suppressed decay 
7r+7r° and doubly-suppressed decay K^tt^ observed experimentally imply a 

sizeable SU(3)-breaking effect. This effect can be explained in the factorization approach 
only if f^^{0) > f^^{0) or R > 1. We find that explanation of the observed ratio R remains 
an unsolved issue in the light-front quark model. 

Not shown in Fig. 5 is the heavy-to-light form factor which is expected to satisfy 

the heavy-quark-symmetry relation at near zero recoil [^]: (/+ + f )^^ ~ l/^m^ and 
(/+ + f-)^'^ ~ l/^/mo. Our light-front calculation shows that in general f^{q^) ~ —f+i.'f) 
is a good approximation for B{D) —>■ it transitions even when is not close to q'^ax' but it 
is only a rough approximation for B{D) K transitions. 

C P Form Factors 

The dependence of the form factors V{q^), ^0,1,2(9^) for B ^ D* transition is depicted 
in Fig. 6. We see that the valence-quark contribution to V, Aq and Ai depends on the choice 
of the or "— " reference frame, corresponding to r(g^) = r+(g^) or r(g^) = r_(g^). In 
general, the form factor in the frame is larger than that in the "— " frame, but they 
become identical at zero recoil where ^'+(q'max) — -i'im&y) — M£)*/Mb [see ( |2.31[ - |2.33|) ]. At 
maximum recoil = 0, we find 

V^^^*(0) = 0.78, A^^*(0) = 0.73, Af^*(0) = 0.68, Af^*(0) = 0.61 , (4.7) 

in the frame where r(0) = r_|_(0) = 1, and 

V^^*(0) = 0.62, <^*(0)=0.58, Af^*(0) = 0.59, ^* (0) = 0.61 , (4.8) 



which allows one to perform an analytic continuation to time-like momentum transfer. The Landau 
singularity there corresponds to our valence-quark contribution, while the non-Landau singularity 
to the non-valence configuration. However, the contribution of the Landau singularity in this 
approach vanishes at the "quark zero recoil" point (see Fig. 14 of 
phenomenon not seen in our direct light-front calculations. 



11] for D ^ K transition), a 
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in the "-" frame where r(0) = r_(0) = {MD*/MBf [see (|]3|)]. As discussed in Sec. II.D, 
no firm predictions for V, Aq, Ai can be made unless the Z-graph contributions are included 
so that they are independent of the "±" frames. Although we do not have a reliable estimate 
for the Z-graph contribution, we know that it is more important for the r = r_ curve than 
the r = Tj^ one. This is because form factors at = do not receive the non-valence 
contribution in the "+" frame because r+(0) = 1. Therefore, (|4.71 ) gives the complete results 
for B ^ D* form factors at = 0. Consequently, the difference between (|4.71) and (4.8) 
must be equal to the non-valence contribution in the "— " frame, namely, 

\/^^*(0) = 0.16, i^^*(0) = 0.15, lf^*(0) = 0.08, lf^*(0)=0. (4.9) 

This implies that for heavy-to-heavy transition, form factors calculated from the valence- 
quark configuration alone and evaluated in the frame should he reliable in a broad kine- 
matic region and become most trustworthy in the close vicinity of maximum recoil. A generic 
feature of the Z-graph effect is illustrated in Fig. 7 by considering the form factor A^^* . As- 
suming thaX the full Aq^* has a dipole behavior shown in Fig. 7 with a pole mass Mpoie = 6.73 
GeV (dash-dotted curve), the difference between the "full curve" and the valence contribu- 
tion should give the non-valence contribution. It is clear that the Z-graph effect in the "— " 
frame (dashed curve) is sizeable in the entire kinemetic region, whereas it is important in 
the frame (solid curve) only when is close to the zero- recoil point. 

For a broad range of g^, we find that Aq^ , v4f , V^^ can be fitted to a dipole form and 
A^^* to a monopole form, in accord with the HQS relations given in ( ^.241) . Experimentally, 
two form-factor ratios defined by 



1- ^ 



V 



BD* 



Ar'{q') 

A^^\q^y 



(4.10) 



{Mb + Md*)\ 

have been extracted by CLEO |^ from an analysis of angular distribution m B ^ D*iu 
decays with the results: 

/2^(g^^^) = 1.18 ±0.30 ±0.12, i?2(gLx) = 0.71 ±0.22 ±0.07. (4.11) 

( p.24| ) implies that, irrespective of the values of g^, -Ri(g^) = -^2(0'^) = 1 in the heavy- 
quark limit. Our light-front calculations yield V^^* (ql^^^) = 1.14, ^f^*(gmax) = 0.83, 
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and ^f'^*(q'max) = 0-96, hence -Ri((?max) = l-H -R2(5'max) = 0.92, in agreement with 
experiment. The predictions of HQET are similar 0: i?i ~ 1.3 ± 0.1 and R2 — 0.8 ± 0.2 . 

As shown in Figs. 8-11, we have also computed the dependence of the form factors for 
B K*, B ^ p, D ^ K* and D ^ p decays. The numerical results of the form factors at 
= are (in the "+" frame): 

B^K*: A^^*(0) = 0.32, Af^*(0) = 0.26 , ^* (0) = 0.23 , y^^*(0) = 0.35 , 
D^K*: A^^*(0) = 0.71, Af^*(0) =0.62, ^*(0) = 0.46 , V^^'(O) = 0.87 , 

B^p: ^^"(0) = 0.28 , Af''(0) = 0.20, Af^(0) = 0.18, \/^''(0) = 0.30 , (4.12) 

D^p: A^^(O) =0.63, ^^''(O) = 0.51 , Af^{0) = OM, l^^^(O) = 0.78 . 

Experimentally, only D K* form factors have been measured with the results [p3| 

V^^* (0) = 1.1 ± 0.2 , ^* (0) = 0.56 ± 0.04 , A^^* (0) = 0.40 ± 0.08 , (4.13) 

obtained by assuming a pole behavior for the dependence. Our predictions for the D K* 
form factors are consistent with experiment. 
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Table II. Form factors for B ^ p and D ^ p transitions at = in various models. 





Reference 


AfP{0) A^P{0) V^P{0) 




Lattice 


BES !§ 
LMMS [|9| 
ELC |0| 
APE P 
UKQCD H 
GSS |3| 


0.22 ± 0.05 0.49 ± 0.21 ± 0.05 0.37 ± 0.11 
0.24 ± 0.12 0.27 ± 0.80 0.53 ±0.31 
0.27tl-tl 0.28t^tt 

nif;+4+22 79+35+10 nP;i+23+9 
"•l"-4-16 L)./Z_35_y U.01_23-6 


0.65+}^+23 0.59;i^}+25 1.07 ±0.49 ±0.35 
0.45 ± 0.04 0.02 ± 0.26 0.78 ± 0.12 

n f^Q+'i'+s n s-^+^o+is 10-1+25+18 


Sum 
Rule 


s 

Ball H 

ABs m 

Narison 
YH |38| 


0.96±0.15 1.21±0.18 1.27±0.12 
0.5 ±0.1 0.4 ±0.2 0.6 ±0.2 
0.24 ±0.04 - 0.28 ±0.06 
0.38 ±0.04 0.45 ± 0.05 0.45 ± 0.05 
0.07 ±0.01 0.16 ±0.01 0.19 ±0.01 


0.5±0.2 0.4±0.1 1.0±0.2 
0.34 ± 0.08 0.57 ± 0.08 0.98 ±0.11 


QM 


ISGW [|9| 
BSW 
Stech @ 
PGM 111 
IV [H 


0.05 0.02 0.27 
0.28 0.28 0.33 
0.32 0.35 0.37 
0.26 ± 0.03 0.31 ± 0.03 0.29 ±0.03 
0.50 0.51 0.70 


0.78 0.92 1.23 

- - 

0.55 0.45 1.08 


LFQM 


this work 

OXT § 
Jaus 
Melikhov (l|] 


0.20 0.18 0.30 

0.21 0.18 0.32 
0.26 0.24 0.35 
0.17-0.26 0.16-0.24 0.22-0.34 


0.51 0.34 0.78 
0.58 0.42 0.93 


HQET 
+ChPT 


CDDGFN [H 
CDDGFN @ 


0.21 0.20 1.04 
0.28 0.19 0.50 


0.55 0.28 1.01 



Form factors for B ^ p and D ^ p transitions at = predicted in various approaches 
(lattice simulations, QCD sum rule, quark model, light-front quark model, and heavy quark 
effective theory together with chiral perturbation theory) are summarized in Table II. We 
have to await further experimental studies in order to test various models. The B ^ K and 
B K* transitions arise from flavor-changing neutral currents induced by QCD corrections. 
It has been found recently that there are two experimental data for B J/ipK^*^ which 
cannot be accounted for simultaneously by all commonly used models Hence, it is 

important to have a reliable estimate of the B K^*^ form factors at = rnj/^ in order to 
test the validity of the factorization approach. Our calculation gives 

Ff^'^im]^^) = 0.66 , Vf^\m]^^) = 0.42 , 

A^'^'im'j/^) = 0.63 , Af^'K/^) = 0.37, = 0.43 , (4.14) 

from valence-quark configuration. 
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As for the dependence of heavy-to-light from factors, we see from Figs. 8-11 that, 
except for V^^ , they all increase with g^, though Ai is flatter than Aq, A2, and V . As we 
have argued before, the valence-quark contribution evaluated in the frame should be 
reliable when is close to maximum recoil. For small g^, we have a dipole behavior for 
Aq, A2, V (except for V^^' and V^^*) and a monopole behavior for Ai, that is, Aq, A2, 
and V increase with faster than Ai. The form factor V for B —>■ p and B — > K* in the 
frame do not have a dipole behavior at small g^ mainly because of the large destructure 
contributions from the Qvk±/Wv terms in (|2.78 ). As a result, the form factor V in B —>■ p 



and B K* decays evaluated in the frame is smaller than that in the "— " frame. 
The g^ dependence of the P ^ V form factors have also been studied in the QCD-sum-rule 
approach with some contradicting results. For example, while is found to decrease from 
g^ = to g^ = 15 GeV^ in |^ (see also P3| , |T7| , |i^ ), such a phenomenon is not seen in PB| , |55| 
(see also Sec. 5.3 of [0]). The sum-rule results of show that the form factors Aq, A2, V 



all have a dipole form while Ai has a monopole form, in accord with our observation. The 
same conclusion is also reached in based on the scaling behavior of heavy-to-light from 
factors in the mg —>■ 00 limit. A recent lattice study of the axial form factors Aq'', Af'' and 
A^'^ |3^ is consistent with the g^ behavior we have in the light-front quark model. 



V. SUMMARY 

The heavy-to-heavy and heavy-to-light form factors in P — P and P ^ V transitions 
are studied in the present paper. In the light-front relativistic quark model, the decay form 
factors are evaluated in a frame where g"*" > and q± = 0, so that it covers the entire 
physical range of momentum transfer and no extrapolation assumption from g^ = or from 
g^ = g^g^^ is required. In previous works using g+ = 0, one can only calculate form factors 
at g^ =0; moreover, the form factors /-(g^) in P — P decay and a-(g^) in P — decay 
cannot be studied. For the first time, we have calculated the P ^ V form factors directly 
at time-like momentum transfers. The main results of this paper are : 

1). We have investigated the behavior of heavy-to-heavy form factors in the heavy-quark 
limit and found that the requirements of heavy-quark symmetry (|2.29| ) for P ^ P transition 



and (p.51| ) for P —>■ V transition are indeed fulfilled by the light-front quark model provided 



that the universal function ({v-v') obtained from P ^ V decay is identical to the Isgur-Wise 
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function ^{v • v') in P — > P decay. 

2) . Contrary to the Isgur-Wise function in P — > P decay, the normahzation of ({v ■ v') 
at zero recoil depends on the hght-front wave function used. We found that the BSW 
amphtude correctly gives .^(1) = 1, but = 0.87. Therefore, this type of wave functions 
cannot describe P ^ V decays in a manner consistent with heavy- quark symmetry. 

3) . Using the Gaussian- type amplitude, the Isgur-Wise function ({v ■ v') has a correct 
normalization at zero recoil and is identical io $,{v ■ v') numerically up to six digits. It can 
be fitted very well with a dipole dependence with Mpoie = 6.65 GeV for P — > D transition. 
However, the predicted slope parameter — 1.24 is probably too large. This may be ascribed 
to the fact that the Gaussian-type wave function does not have enough high-momentum 
components at large k^. 

4) . The valence-quark and non- valence contributions to form factors are in general depen- 
dent on the recoiling direction of the daughter meson relative to the parent meson, but their 
sum should not. Although we do not have a reliable estimate of the pair-creation effect, we 
have argued that, for heavy-to-heavy transition, form factors calculated from the valence- 
quark configuration evaluated in the frame should be reliable in a broad kinematic 
region, and they become most trustworthy in the vicinity of maximum recoil. 

5) . The form factors Fi, Aq, A2, V (except for V^^ and V^^*) all exhibit a dipole 
behavior, which Fq and Ai show a monopole behavior in the close vicinity of maximum 
recoil for heavy-to-light transition, and in a broader kinematic region for heavy-to-heavy 
decays. Therefore, Pi, Aq,A2, V increase with faster than Pq and A^. 
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FIGURE CAPTIONS 



Fig. 1 The Feynman triangle diagram and the corresponding hght-front subdiagrams. 
Diagram (a) corresponds to the valence-quark configuration and diagram (b) to the non- 
valence configuration. 

Fig. 2 The Isgur-Wise function ^ {v ■ v') for P ^ P transition calculated using Gaussian- 
type (sohd fine) and BSW-type (dashed line) hght-front wave functions. For comparison, a 
curve ior 1/v ■ v' is also shown. 

Fig. 3 The Isgur-Wise function C(^"^') for P — transition calculated using Gaussian- type 
(solid line) and BSW-type (dashed line) light-front wave functions. 

Fig. 4 Form factors Pi(g^) = f+{(f) and Po(q'^) for P ^ D transition arising from the 
valence-quark configuration. Dashed curves are fits to Pi in a dipole form with Mpoie = 6.59 
GeV and to Pq in a monopole form with Mpde = 7.90 GeV. 

Fig. 5 The form factor /+(g^) for P — > tt, P — > P — > tt and P — > X transitions arising 
from the valence-quark configuration. 

Fig. 6 Form factors V{q^), AQ{q^), Ai{q^) and A2{q^) for B ^ D* transition. Solid fines 
are the valence contribution evaluated in the frame where r(g^) = r+(g^), and dashed 
lines in the "— " frame where r(g^) = r_(g^). The contribution to the form factor A2 is 
independent of the choice of the or "— " frame. 

Fig. 7 An illustration of the non- valence contribution to the form factor A^^* , whose 
general feature applies to Af^* and V^^* as well. Valence-quark contributions to Aq^' 
evaluated in the frame (sohd fine) and in the "— " frame (dashed fine) are the same 
as in Fig. 6. The corresponding non-valence contribution are extracted in respective frames 
by assuming that the full Aq'^* has a dipole behavior with Mpoie = 6.73 GeV (dash-dotted 
line) . 

Fig. 8 Same as Fig. 6 except for P — > transition. 
Fig. 9 Same as Fig. 6 except for P — > p transition. 
Fig. 10 Same as Fig. 6 except for D — > K* transition. 
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Fig. 11 Same as Fig. 6 except for — > p transition. 
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